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Abstract 

We suggest an effective approach to separation of variables in the Schrodinger 
equation with two space variables. Using it we classify inequivalent potentials 
V(xi,X2) such that the corresponding Schrodinger equations admit separation of 
variables. Besides that, we carry out separation of variables in the Schrodinger 
equation with the anisotropic harmonic oscillator potential V = k\x\ + k^x\ and 
obtain a complete list of coordinate systems providing its separability. Most of 
these coordinate systems depend essentially on the form of the potential and do not 
provide separation of variables in the free Schrodinger equation (V = 0). 



I. Introduction 

The problem of separation of variables [SV) in the two-dimensional Schrodinger equation 

iu t + u XlXl +u X2X2 = V(x 1 ,x 2 )u (1) 

as well as the most of classical problems of mathematical physics can be formulated in 
a very simple way (but this simplicity does not, of course, imply an existence of easy 

*On leave from the Institute of Mathematics of the Academy of Sciences of Ukraine, Tereshchenkivska 
Str.3, 252004 Kiev, Ukraine 

e-mail: asrz@pta3.pt.tu-clausthal.de 



way to its solution). To separate variables in Eq.® one has to construct such functions 
R(t,x), u>i(t,x), u 2 (t,x) that the Schrodinger equation (P after being rewritten in the 
new variables 

Z = t, Z 1 =L0 1 (t,x), Z 2 =U 2 (t,x), ^ 

v(z , z) = R(t, x)u(t, x) 

separates into three ordinary differential equations (ODEs). From this point of view the 
problem of SV in Eq.(^) is studied in 

But no less important problem is the one of description of potentials V(xi,x 2 ) such 
that the Schrodinger equation admits variable separation. That is why saying about SV 
in Eq. (HD we imply two mutually connected problems. The first one is to describe all such 
functions V(x\, x 2 ) that the corresponding Schrodinger equation ([]]) can be separated into 
three ODEs in some coordinate system of the form (|2]) (classification problem). The second 
problem is to construct for each function V(xi,x 2 ) obtained in this way all coordinate 
systems ® enabling us to carry out SVin Eq. ([[]). 

Up to our knowledge, the second problem has been solved provided V = |fj, || and 

V 2 + P%2 2 E- The first 

one was considered in a restricted sense in Iffl. Authors 
using symmetry approach to classification problem obtained some potentials providing 
separability of Eq. QTj) and carried out .SVin the corresponding Schrodinger equation. But 
their results are far from being complete and systematic. The necessary and sufficient 
conditions imposed on the potential V(xi,x 2 ) by the requirement that the Schrodinger 
equation admits symmetry operators of an arbitrary order are obtained in ||. But so far 
there is no systematic and exhaustive description of potentials V(xi,x 2 ) providing SVin 
Eq.©. 

To have a right to say about description of all potentials and all coordinate systems 
making it possible to separate the Schrodinger equation one has to give a definition of SV. 
One of the possible definitions of SVin partial differential equations (PDEs) is proposed in 
our paper || . It is based on the concept of Ansatz suggested by Fushchych |7j and on ideas 
contained in the paper by Koornwinder 0. The said definition is quite algorithmic in 
the sense that it contains a regular algorithm of variable separation in partial differential 
equations which can be easily adapted to handle both linear || |9| and nonlinear [fLOf 



PDEs. In the present paper we apply the said algorithm to solve the problem of SV in 
Eq.®. 

Consider the following system of ODEs: 

Ui{uju (pi, -: — ;Ai,A 2 ), (3) 



d 2 ip 2 dip 



U 2 (uj 2 ,p 2 ,- — ;Ai,A 2 ), 



duo dto 2 
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where Uq, U\, U2 are some smooth functions of the corresponding arguments, Ai, A2 C R 1 
are arbitrary parameters (separation constants) and what is more 



rank 



= 2 (4) 

fi=0 a=l 



d\ a 

(the last condition ensures essential dependence of the corresponding solution with sepa- 
rated variables on Ai, A2, see ||). 

Definition 1. We say that Eq.(|l|) admits 5" Fin the system of coordinates t, ui(t,x), 
dj 2 (t, x) if substitution of the Ansatz 

u = Q(t, xjipofyipx^uxit, x))tp 2 (u2(t, x)) (5) 

into (|l|) with subsequent exclusion of the derivatives ^r, *^jr according to equations 

(|) yields an identity with respect to <p , <p lt ip 2 , Ai, A 2 . 

Thus, according to the above definition to separate variables in Eq.(|l|) one has 

• to substitute the expression (|5|) into ([[]), 

• to exclude derivatives ^r, with the help of equations (H), 

• to split the obtained equality with respect to the variables (p , (fx, <p 2 , Ai, 
A2 considered as independent. 

As a result one gets some over-determined system of PDEs for the functions Q(t,x), 
u>i(t,x), uj 2 (t,x). On solving it one obtains a complete description of all coordinate 
systems and potentials providing SVin the Schrodinger equation. Naturally, an expression 
complete description makes sense only within the framework of our definition. So if one 
uses a more general definition it may be possible to construct new coordinate systems and 
potentials providing separability of Eq. ([]]). But all solutions of the Schrodinger equation 
with separated variables known to us fit into the scheme suggested by us and can be 
obtained in the above described way. 



II. Classification of potentials V(x\,X2). 

We do not adduce in full detail computations needed because they are very cumbersome. 
We shall restrict ourselves to pointing out main steps of the realization of the above 
suggested algorithm. 

First of all we make a remark, which makes life a little bit easier. It is readily seen 
that a substitution of the form 



Q -> Q' = Q*xM*2M, 

u a -> uj' a = n o (o; ), a = 1,2, (6) 
A a -> A^ = A a (Ai, A 2 ), a = 1,2, 
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does not alter the structure of relations (J3|), (£|), (||). That is why, we can introduce the 
following equivalence relation: 



(wi, u; 2 , Q) ~ u)' 2 , Q') 

provided (|j) holds with some \l/ a , f2 a , A a . 

Substituting @ into (|l|) and excluding the derivatives with the use of 

equations ([D we get 

i(Qt¥>o<Pi<P2 + QU tpiip2 + QuuVq<PW2 + QuzttpoVifa) + (AQ)ip ^1^2 
+2Q Xa uj 1Xa ip ip 1 (p 2 + 2Qx a UJ 2x a ¥o<Pi ( p2 + C?((Aa;i)^o0i^2 + (Au; 2 )<^ <£W2 

where the summation over the repeated index a from 1 to 2 is understood. Hereafter an 
overdot means differentiation with respect to a corresponding argument and A = d% +d% . 

Splitting the equality obtained with respect to independent variables (fx, y2 2 , , ffi, 
Ai, A 2 we conclude that ODEs (H) are linear and up to the equivalence relation (]6|) can 
be written in the form 

(\ 1 R 1 (t)+\ 2 R 2 (t) + Ro(t))(p , 
(XiBixfa) + \ 2 B 1 2(uj 1 ) + Box{ui))ipu 

\\xB 2 x{u)2) + A 2 522(W 2 ) + B 02 {UJ2))^2 

oji, &2, Q satisfy an over-determined system of nonlinear 



1- ^lx b ^2x b — 0, 

2. B la (LOx)L0 1Xb uJxx b + £20(^2)^2^2*6 + Ra{t) = 0, a = 1, 2, 

3. 2^^^ + Q(«u; a t + Aw a ), a = 1, 2, (7) 

4. (5 i(^i Viz b ^Lr b + £02(^1)^2x^2x6) Q + iQt + AQ + Ro{t)Q 
-V(xx,x 2 )Q = 0. 

Thus, to solve the problem of S'Ffor the linear Schrodinger equation it is necessary to 
construct general solution of system of nonlinear PDEs (0). Roughly speaking, to solve a 
linear equation one has to solve a system of nonlinear equations] This is the reason why 



d(fo 

i>—r- 
fit 

<P<Pi 

dwf 

d?<P2 

dio\ 

and what is more, functions 
PDEs 
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so far there is no complete description of all coordinate systems providing separability of 
the four- dimensional wave equation ||. 

But in the case involved we have succeeded in integrating system of nonlinear PDEs 
(0). Our approach to integration of it is based on the following change of variables 
(hodograph transformation) 



z = t, zi = Z 1 {t,u u u 2 ), z 2 = Z 2 {t,oj u u> 2 ), 

V\ = Xt, v 2 = x 2 

where zq, Zi, z 2 are new independent and v±, v 2 are new dependent variables correspond- 
ingly. 

Using the hodograph transformation determined above we have constructed the gen- 
eral solution of equations 1-3 from (0). It is given up to the equivalence relation (|) by 
one of the following formulae: 



wi = A(t) Xl + Wi(t), u 2 = B{t)x 2 + W 2 {t), 

( i ( A 2 B 2 x i f Wx W 2 



Q(t, x) = exp|-- {— Xl + —x 2 ) - - + — x 2 

2. u)\ — - \n(x\ + I?) + W(t), u> 2 = arctan — , 

2 x 2 

Q(t,x) = exp|-^(a^ + x|)|; 

3. xi = ^(t)K 2 -cj 2 2 ) + ^i(t), x 2 = W(t)u 1 u 2 + W 2 (t), (8) 

Q(t,x) = exp|^((a; 1 - Wrf + (x 2 - W 2 f) + % -{W x x x + fe 2 )|; 

4. xi = W(t) cosh ui cos uj 2 + Wi(t), x 2 = W(t) sinho;i sino; 2 + W 2 (t), 
Q(t,x) = exp|^((a; 1 - Wrf + (x 2 - W 2 f) + % -{W lXl + fe 2 )|; 

Here A, B, W, W±, W 2 are arbitrary smooth functions on t. 

Substituting the obtained expressions for the functions Q, u>i, lo 2 into the last equa- 
tion from the system (|7|) and splitting with respect to variables xi, x 2 we get explicit 
forms of potentials V(xi,x 2 ) and systems of nonlinear ODEs for unknown functions 
A(t), B{t), W(t), Wi(t), W 2 (t). We have succeeded in integrating these and in con- 
structing all coordinate systems providing SVin the initial equation ([!]). 

Here we consider in detail integration of the fourth equation of system (^) for the case 
2 from (H), since computations needed are not so lengthy as for other cases. 

First, we make several important remarks which introduce an equivalence relation on 
the set of potentials V(xx,x 2 ). 

Remark 1. The Schrodinger equation with the potential 
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V(x±, x 2 ) = hxi + k 2 x 2 + k 3 + Vi{k 2 xi - kix 2 ), (9) 

where k\, k 2 , k 3 are constants, is transformed to the Schrodinger equation with the 
potential 

V'(x[,x' 2 ) = V 1 {k 2 x' 1 - k lX ! 2 ) (10) 
by the following change of variables: 

t — it ; x — — x ~f~ t k j 

(i (ii) 

w' = ■uexp|-(A;^ + A;^)^ 3 + it(k%Xi + k 2 x 2 ) + 2A; 3 t|. 

It is readily seen that the class of Ansatze is transformed into itself by the above 
change of variables. That is why, potentials and (|10|) are considered as equivalent. 
Remark 2. The Schrodinger equation with the potential 



V( Xl ,x 2 ) = k{x\ + xl) + ^1(^)^1 + 



'12) 



,x 2 . 

with k = const is reduced to the Schrodinger equation with the potential 



V'(x 1 ,x 2 )=V 1 (^)(x' 1 2 + x' 1 2 )- 1 (13) 



X' 



by the change of variables 

t' = a(t), a? — p(t)x, u = uexp^(t)(xl + x\) + 5(t)}, 
where (a(t), (3(t), 7(t), S(tfj is an arbitrary solution of the system of ODEs 

7 - 4 7 2 = k, $ - 47/3 = 0, a - P 2 = 0, 5 + 4 7 = 

such that [3 7^ 0. 

Since the above change of variables does not alter the structure of the Ansatz (|5|), when 
classifying potentials V(xi,x 2 ) providing separability of the corresponding Schrodinger 
equation we consider potentials dl2[) , flilf ) as equivalent. 

Remark 3. It is well-known (see e.g. [ |TT| , |12|) that the general form of the invariance 
group admitted by Eq.(|l|) is as follows 

t' = F{t,8), x' a = g a {t,x,9), a = 1,2, 
v! = h(t, x, 6)u + U (t, x) , 

where 6 = (#1, 6 2 , . . . , 9 n ) are group parameters and U{t,x) is an arbitrary solution of 
Eq.®. 
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The above transformations also do not alter the structure of the Ansatz ([5]). That is 
why, systems of coordinates f, x[, x 2 and t, x±, x^ are considered as equivalent. 

Now we turn to the integration of the fourth equation of system ([?]). Substituting into 
it the expressions for the functions Ui, uj 2 , Q given by formulae 2 from (^j we get 

V( Xl ,x 2 ) = (B 01 (u 1 ) + B 02 (u; 2 ))exp{-2(u; 1 - W)} + -(W -W 



4 V ' (14) 

x exp{2(w! -W)} + Ro{i) - iW. 

In the above equality £> i, B 02 , Ro(t), W(t) are unknown functions to be determined 
from the requirement that the right-hand side of flT4] ) does not depend on t. 
Differentiating ( |i~4"|) with respect to t and taking into account the equalities 

u lt = W, u 2t = 

we have 

VFexp{-2(^! - W)}B 01 + a(t) exp{2(cj! - W)} + (3(t) = 0, (15) 

where a(t) = \{W - W 2 ), (3{t) = R - iW . 

Cases W = and W ^ have to be considered separately. 
The case 1. W = 0. 

In this case W = C = const, Ro = 0. Since coordinate systems u>i, u 2 and u>i + 
Ci, UJ2 + Ci are equivalent with arbitrary constants Ci, C2, choosing G\ = — C, C2 = 
we can put C = 0. Hence it immediately follows that 



V{xi,x 2 ) 



?oi Q In (a^ + 0^2) j + -B02 ^arctan ^ 

where i?oi, -B02 are arbitrary functions. And what is more, the Schrodinger equation ([!]) 
with such potential separates only in one coordinate system 



lui = -m(x 1 + x 2 ), co>2 = arctan — . (16) 
2 £2 

TTie case 2. W ^ 

Dividing ( [14] ) into iyexp{— 2(co>i — W)} and differentiating the equality obtained with 
respect to t we get 

exp{4a; 1 }^-(d(M/)- 1 exp{-4H/}) + exp{2u 1 }^/3(W)- 1 exp{-2W}) =0, 

tLL Lib 

whence 

— (d(W)- 1 exp{-4W}) = 0, -^-(/5(H/)- 1 exp{-2H/}) = 0. 
Integration of the above ODEs yields the following result: 
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a = d exp{4W} + C 2 , P = C 3 exp{2W} + C 4 , 

where C 3 , j — 1,4 are arbitrary real constants. 

Inserting the result obtained into (|THj) we get an equation for £? 01 

JBqi = -4Ciexp{4u;i} - 2C 3 exp{2u;i}, 
which general solution reads 

Bqi = —Ci exp{4a;i} - C 3 exp{2a;i} + C 5 . 
In the above equality C5 is an arbitrary real constant. 

Substituting the expressions for a, (3, B m into Eq.flUP we have the explicit form of 
the potential V(xi,x 2 ) 



V(x 1 ,x 2 ) 



Bq2 [ arctan — ) + C 5 



x{ + x*)- 1 + C 2 (xj + x 2 2 ) + C 4; 



x 2 . 

where Bq 2 is an arbitrary function. 

By force of the Remarks 1, 2 we can choose C 2 = C4 = 0. Furthermore, due to 
arbitrariness of the function B 02 we can put C5 = 0. 

Thus, the case W 7^ leads to the following potential: 

V(xi,x 2 ) = B 02 (arctan — ) (x\ + xl)^ 1 . (17) 
V x 2 J 

Substitution of the above expression into Eq. flliD yields second-order nonlinear ODE 
for the function W = W(t) 

W-W 2 = Ad exp{AW}, (18) 
while the function Rq is given by the formula 

R = iW + C 3 exp{2W}. 

Integration of ODE ( p!8|) is considered in detail in the Appendix 1. Its general solution 
has the form 
under G\ 7^ 

W = -- In ((at - bf - 4Ci) + \ In a, 

under C\ = 

W = a-\n(t + b). 

Substituting obtained expressions for W into formulae 2 from (§) and taking into 
account the Remark 3 we arrive at the conclusion that the Schrodinger equation ([!]) 
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with the potential ( |TTD admits SV in two coordinate systems. One of them is the polar 
coordinate system ( |i~6|) and another one is the following: 

11 X\ 

Ui = - In (x\ + x\) In [t 2 ±1), uj 2 = arctan — . (19) 

2 2 X2 

Consequently, the case 2 from (|8|) gives rise to two classes of the separable Schrodinger 
equations ([I]). 

Cases 1, 3, 4 from (|8|) are considered in an analogous way but computations involved 
are much more cumbersome. As a result, we obtain the following list of inequivalent 
potentials V(xi,x 2 ) providing separability of the Schrodinger equation. 

1. V(x u x 2 ) = V 1 (x l ) + V 2 (x 2 ); 

(a) V(xi, x 2 ) = k x x\ + k 2 Xi 2 + V 2 (x 2 ), k 2 ^ 0; 

i. V(xi,x 2 ) = k\x\ + k 2 x\ + k 3 Xi 2 + k^x^ 2 , k^k^ ^ 0, 
k\ + k\± 0, h ^ k 2 - 

ii. V(xi, x 2 ) = k x x\ + k 2 x± 2 , k x k 2 ^ 0; 
hi. V(x l , x 2 ) = k x x\ 2 + k 2 x 2 2 ; 

(b) V(xi,x 2 ) = k x x\ + V 2 (x 2 ); 

i. V(xi, x 2 ) = kix\ + k 2 x\ + /c 3 x2 2 , fcifc 3 7^ 0, fci/ fc 2 ; 

ii. V(x 1 ,x 2 ) = k x x\ + k 2 x\, k x k 2 ^ 0, fci / fc 2 ; 
hi. y(xi,X2) = k\x\ + k 2 x 2 2 , k\ 7^ 0; 

2. V(xi, i 2 ) = V x {x\ + x|) + F 2 (—) (xl + xl)' 1 ; 

\x 2 J 

(a) V(x 1 ,x 2 ) = V 2 (^yx 2 1 + x 2 2 )- 1 ; 

(b) V(xi, x 2 ) = h(x 2 + x 2 )- 1 ' 2 , h ± 0; 

3. V{x 1 ,x 2 ) = [v 1 {u 1 ) + V 2 {u 2 )){lo 2 + lo 2 )-\ 
where uj\ — uj 2 = 2x\, uj\uj 2 = x 2 ; 

4. V(x 1 ,x 2 ) = (VxM + F 2 (^2))(sinh 2 u) x + sin 2 uj 2 )~\ 
where coshu;! coscj 2 = X\, sinho>i sinu^ = x 2 ; 

5. V(x ll x 2 ) = 0. 
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In the above formulae Vi, V2 are arbitrary smooth functions, A4, k 2 , k 3 , k 4 are arbi- 
trary constants. 

It should be emphasized that the above potentials are not inequivalent in a usual sense. 
These potentials differ from each other by the fact that the coordinate systems providing 
separability of the corresponding Schrodinger equations are different. As an illustration, 
we give the following figure: 




Fig.l 
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where r = (x\ + x 2 Y^ 2 and by the symbol V^\j = 1,4 we denote the potential given in 
the above list under the number j. Downarrows in the Fig.l indicate specifications of the 
potential Vxi,x 2 ) providing new possibilities to separate the corresponding Schrodinger 
equation (JJ). 

The Schrodinger equation ([!]) with arbitrary function V(xi,x 2 ) (level 1 of the Fig.l) 
admits no separation of variables. Next, Eq. (Jl|) with the "root" potentials V^' (level 2), 
Vi, V 2 being arbitrary smooth functions, separates in the Cartesian (j = 1), polar (j = 2), 
parabolic (j = 3) and elliptic (j = 4) coordinate systems, correspondingly. Specifying the 
functions Vi, V 2 (i.e. going down to the lower levels) new possibilities to separate variables 
in the Schrodinger equation ([!]) arise. For example, Eq.([I]) with the potential V 2 (^)r~ 2 , 
which is a particular case of the potential V^ 2 \ separates not only in the polar coordinate 
system (|I1J) but also in the coordinate systems (|i9|). The Schrodinger equation with the 
Coulomb potential fcir -1 , which is a particular case of the potentials V^ 2 \ separates 
in two coordinate systems (namely, in the polar and parabolic coordinate systems, see 
below the Theorem 4). An another characteristic example is a transition from the po- 
tential V^ 1 ) to the potential k%x\ + V 2 (x 2 ). The Schrodinger equation with the potential 
admits SV in the Cartesian coordinate system ojq = t, uj\ — xi, co 2 = x 2 only, while 
the one with the potential k\x\ + V 2 (x 2 ) separates in seven {ki < 0) or in three (ki > 0) 
coordinate systems. 

A complete list of coordinate systems providing 5" Fin the Schrodinger equations with 
the above given potentials takes two dozens of pages. Therefore, we restrict ourself 
to considering the Schrodinger equation with anisotropic harmonic oscillator potential 
V(xi, x 2 ) = kixf + k 2 x\, k\ 7^ k 2 and Coulomb potential V(x\, x 2 ) = k\(x\ + x 2 )^ 1 / 2 . 

III. Separation of variables in the Schrodinger equa- 
tion with the anisotropic harmonic oscillator and the 
Coulomb potentials. 

Here we will obtain all coordinate systems providing separability of the Schrodinger equa- 
tion with the potential V(xi,x 2 ) = k\x\ + k 2 x\ 

iu t + u XlXl + u X2X2 = (kix\ + k 2 x 2 2 )u. (20) 

In the following, we consider the case ki ^ k 2 , because otherwise Eq.(|l|) is reduced 
to the free Schrodinger equation (see the Remark 2) which has been studied in detail in 
[1-3]. 

Explicit forms of the coordinate systems to be found depend essentially on the signs 
of the parameters ki, k 2 . We consider in detail the case, when ki < 0, k 2 > (the cases 
k\ > 0, k 2 > and ki < 0, k 2 < are handled in an analogous way). It means that 
Eq.(P0"D can be written in the form 
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iu t + u XlXl + u X2X2 + -{a 2 x\ - b 2 x\)u = 0, (21) 

where a, b are arbitrary non-null real constants (the factor \ is introduced for further 
convenience) . 

As stated above to describe all coordinate systems t, cux(t,x), u) 2 (t, x) providing 
separability of Eq.(^) one has to construct the general solution of system (H) with 
V(xi,x 2 ) = — \{a 2 x\ — b 2 x 2 ). The general solution of equations 1-3 from (|7|) splits into 
four inequivalent classes listed in (H). Analysis shows that only solutions belonging to the 
first class can satisfy an equation 4 from (^). 

Substituting the expressions for u x , uj 2 , Q given by the formulae 1 from (§) into the 
equation 4 from (0) with V(x x , x 2 ) = — \{a 2 x\ — b 2 x 2 ) and splitting with respect to x±, x 2 
one gets 



-Boi(^i) = oiiu{ + a 2 u u £02(^2) = fixu\ + f3 2 u 2 , 



A! - \A 



2 



4 ai A 4 + a 2 = 0, (22) 



B \ -(^) 2 -^-b 2 = 0, (23) 



sy \BJ 

9 X - 26 x - - 2(2a 1 6 1 + a 2 )A A = 0, (24) 
#2 -2^2^- 2(2/3A + /? 2 ) J B 4 = 0. (25) 
Here «i, 02, are arbitrary real constants. 

Integration of the system of nonlinear ODEs ( ^2j) -fl25|) is carried out in the Appendix 
1. Substitution of the formulae (A. 2), (A.4)-(A.6), (A.8)-(A.ll) into the corresponding 
expressions 1 from (||) yields a complete list of coordinate systems providing separability 
of the Schrodinger equation ([H]). These systems can be transformed to canonical form if 
we use the Remark 3. 

The invariance group of Eq. Q2"ip is generated by the following basis operators ]TT| : 



P =d t , I = ud u , M = md u , Qoo = U(t, x)d u 
%a 

Pi = cosh at d Xl + ~^{ x i sinh at)ud u , 
%b 

P 2 = co$btd X2 — — (x2 sin bt)ud u , (26) 

in 

G\ = sinh at d Xl + ~w( x i cosh at)ud u , 
%b 

G2 = sinbtd X2 + — (x 2 co$bt)ud u , 
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where U(t, x) is an arbitrary solution of Eg. (|2T|) . 

Using the finite transformations generated by the infinitesimal operators ( p6|) and the 
Remark 3 we can choose in the formulae (A.4)-(A.6), (A. 8), (A.10), (A. 11) C 3 = C 4 = 
Di = 0, D 3 = D 4 = 0, C 2 = -D2 = 1- As a result, we come to the following assertion. 

Theorem 1. The Schrddinger equation ( 21 ) admits SV in 21 inequivalent coordinate 
systems of the form 



Uq = t, 



u x (t,x) 



(27) 



where U\ is given by one of the formulae from the first and u 2 - by one of the formulae 
from the second column of the Table 1. 

Table 1. 





uj 2 (t,x) 


a;i(sinha(t + C)) +a(sinh a(t + C)\ 2 


x 2 (smbt)- 1 + (3(smbt)- 2 


xi^cosha(t + C)^j +a(^cosh a(t + C)J 
x\ exp(±ai) + a exp(±4at) 

X\ (a + sinh 2a(t + C)j ^ 
Xi (a + cosh 2a(t + C)\ ^ 
Xi (a + exp(±2at)^) 


x 2 (P + sm2bt)- 1 / 2 
x 2 




Xi 





Here C, a, (3 are arbitrary real constants. 

There is no necessity to consider specially the case when in fl2"oD k% > 0, k 2 < 0, 
since such an equation by the change of independent variables u(t, x\,x 2 ) — * u(t,x 2 ,xi) 
is reduced to Eq.(|2l|). 

Below we adduce without proof the assertions describing coordinate systems providing 
SVm Eq.(g(]) with k x < 0, k 2 < and k x > 0, k 2 > 0. 

Theorem 2. The Schrddinger equation 



lUt + u XlX1 + u 



,,., ;2 + -^{a 2 xl + b 2 x 2 2 )u = 



(28) 



with a 2 7^ 4b 2 admits SV in 49 inequivalent coordinate systems of the form where lo\ 
is given by one of the formulae from the first and u 2 - by one of the formulae from the 
second column of the Table 2. Provided a 2 = 4b 2 one more coordinate system should be 
included into the above list, namely 
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Co> — t, U>1 — Uj\ = 2X\, UJ\UJ 2 = X 2 . 



(29) 



Table 2. 



-i 



x\ (sinh a(t + C) j +ct(sinh a(t + C) 

xi(cosha(£ + C)) +a(cosha(t + C) 
Xi exp(±a£) + aexp(±4at) 

Xi (a + sinh 2a(t + C)) 

xi (a + cosh 2a(£ + C)j 

X\ (a + exp(±2a£)) ^ 

Xi 



cj 2 (*,x) 

^(sinhfrfp 1 + P(smhbty 2 

x 2 (cosh bt)" 1 + /3(coshto)~ 2 
x 2 exp(±6t) + /3exp(±46t) 

x 2 (/? + sinh26t)- 1 / 2 

x 2 (/5 + cosh2to)- 1 / 2 

x 2 (/3 + exp(±2fc) A 
x 2 



Here C, a, (3 are arbitrary constants. 
Theorem 3. The Schrddinger equation 



iu t + u xlXl + u 



X2X2 



1 I 2 2 , t2 2\ 

-(a x l + b x 2 )u 







(30) 



tt>i£/i a 2 46 2 admits SV in 9 inequivalent coordinate systems of the form $1%, inhere uj\ 
is given by one of the formulae from the first and u 2 - by one of the formulae from the 
second column of the Table 3. Provided a 2 = 4f> 2 ; the above list should be supplemented 



by the coordinate system ft29]). 



Table 3. 



UJi(t,x) 


L0 2 (t,x) 


Xi(sina(t + C)) +a(sin a(t + C)) 
xi(^ + sin 2a(t + C))~ V2 

Xi 


x 2 (sin&£) -1 + [3(smbt)~ 2 

x 2 (/3 + sin26t)- 1 / 2 
x 2 
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Here C, a, (3 are arbitrary constants. 

Remark 4- If we consider (1) as an equation for a complex- valued function u of three 
complex variables t, xi, x 2 , then the cases considered in the Theorems 1-3 are equivalent. 
Really, replacing, when necessary, a with la and b by ib we can always reduce Eqs.(pT|), 
(|8|) to the form (Q). It means that coordinate systems presented in the Tables 1, 2 are 
complex equivalent to those listed in the Table 3. But if u is a complex-valued function 
of real variables t, Xi, x 2 it is not the case. 

Theorem 4. The Schrddinger equation with the Coulomb potential 

iu t + u XlXl + u X2X2 - k x {xl + x 2 2 y 1/2 u = 
admits SV in two coordinate systems ^Tdj), $2$). 



It is important to note that explicit forms of coordinate systems providing separability 



of Eqs.(piD, fl28|) , (|30|) depend essentially on the parameters a, b contained in the potential 
V(xi,x 2 ). It means that the free Schrddinger equation (V = 0) does not admit SV in 
such coordinate systems. Consequently, they are essentially new. 



IV. Conclusion 

In the present paper we have studied the case when the Schrddinger equation (|l|) separates 
into one first-order and two second-order ODEs. It is not difficult to prove that there are 
no functions Q(t,x), oJ^t, x), fi = 0,2 such that the Ansatz 

u = Q(t, x)(p (iJo(t, x^ip^u^t, x)^(p 2 (u 2 (t, x)^ 

separates Eq.(|I|) into three second-order ODEs (see Appendix 2). Nevertheless, there 
exists a possibility for Eq.([l|) to be separated into two first-order and one second-order 
ODEs or into three first-order ODEs. This is a probable source of new potentials and new 
coordinate systems providing separability of the Schrddinger equation. It should be said 
that separation of the two-dimensional wave equation 

Utt ~ u xx = V(x)u 

into one first-order and one second-order ODEs gives no new potentials as compared with 
separation of it into two second-order ODEs. But for some already known potentials new 
coordinate system providing separability of the above equation are obtained 0. 

Let us briefly analyze a connection between separability of Eq.(|l]) and its symmetry 
properties. It is well-known that each solution of the free Schrddinger equation with 
separated variables is a common eigenfunction of two mutually commuting second-order 
symmetry operators of the said equation [Q, |3|. And what is more, separation constants 
Ai, A 2 are eigenvalues of these symmetry operators. 

We will establish that the same assertion holds for the Schrddinger equation Let 
us make in Eq.f]]]) the following change of variables: 
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u = Q(t, x)U (t, ui(t, x) , u> 2 (t, xj) , (31) 

where (Q, cui, uo 2 ) is an arbitrary solution of the system of PDEs (|7|). 

Substituting the expression into (|l|) and taking into account equations (|7|) we get 

Q\iU t + (U WlUJl - B Ql (uJi)U^uJi Xa uJi Xa + (U U2W2 - B 02 (u 2 )U^u 2Xa L0 2Xa j = 0. (32) 

Resolving equations 2 from the system (|7|) with respect to uj\ Xa uJi Xa and uj 2xa u 2xa we 
have 

u 2Xa u 2Xa = -(Rx^Bx^ui) - R 2 {t)B n {u) 1 ) 

where 5 = Bu{uji)B 22 {uj 2 ) — Bi 2 (u>i)B 2 i(uj 2 ) (5 7^ by force of the condition 

Substitution of the above equalities into Eq.fl32"l) with subsequent division by Q 7^ 
yields the following PDE: 



iU t + ^f-(B 12 M(U^ 2 - B 02 (u 2 )u) - B^u^U^ - B 01 MU 

+ ^-(b^m^u^ - Soi(wi)^) - - B m {u 2 )u)"j = 0. 



(33) 



Thus, in the new coordinates t, u>i, cu 2 , U(t,u>i,u> 2 ) Eq.(|I]) takes the form (|3~3"D. 
By direct (and very cumbersome) computation one can check that the following 
second-order differential operators: 



* = ^K-B m M)-^(dl 2 -B« 2 (u 2 
X 2 = -^(dl i -B 0lM ) + ^(dl 2 -B Q2M 



commute under arbitrary B 0a , B ab , a, b = 1, 2, i.e. 

X u X 2 ] = X 1 X 2 - X 2 X X = 0. (34) 
After being rewritten in terms of the operators X±, X 2 Eq.([£J) reads 

(id t - Ry[t)Xx - R 2 {t)X 2 )u = 0. 

Since the relations 
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i3 t -R 1 {t)X 1 -R 2 {t)X 2 , X a 



0. 



1,2 



(35) 



hold, operators Xi, X2 are mutually commuting symmetry operators of Eq. (p3|) . Fur- 
thermore, solution of Eq.(|3~3"D with separated variables U = <po(t) ^i(^i) ^2(^2) satisfies 
the identities 



X a U = X a U, a =1,2. 



(36) 



Consequently, if we designate by X[, X' 2 the operators X x , X 2 written in the initial 
variables t, x, u, then we get from (R3j)-(R6]) the following equalities: 



id t + A - V(x 1 ,x 2 ), X' a ] =0, a =1,2 
x[, x' 2 = 0, 

X' a u = \ a u, a = 1,2. 



where u = Q(t,x)ip {t)^i(uJi)(p 2 (uj 2 ). 

It means that each solution with separated variables is a common eigenfunction of 
two mutually commuting symmetry operators X[, X 2 of the Schrodinger equation ([!]), 
separation constants Ai, A 2 being their eigenvalues. 

Detailed study of the said operators as well as analysis of separated ODEs for functions 
fin A* = 0, 2 (in the way as it is done for the free Schrodinger equation in || [|) is in 
progress and will be a topic of our future publications. 
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Appendix 1. Integration of nonlinear ODEs (p2|)— (p5|). 

Evidently, equations (p2|)-(p5|) can be rewritten in the following unified form: 

2 



— Aay = k, 



z - 2z- - 2(2az + (3)y 4 

y 



0. 



(Al) 



Provided k = — a 2 < 0, system (A.l) coincides with equations (p2|) , (p4|) and under 
k = b 2 > - with equations ([23]), @. 

First of all, we note that the function z = z(t) is determined up to addition of an 
arbitrary constant. Really, the coordinate functions u a have the following structure: 
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u a = yx a + z, a = 1,2. 

But the coordinate system t, u>i, u 2 is equivalent to the coordinate system t, lo\ + 
Ci, u 2 + C 2 , C a E R 1 . Hence it follows that the function z(t) is equivalent to the function 
z(t) + C with arbitrary real constant C. Consequently provided a ^ 0, we can choose in 
(A.l) = 0. 

The case 1. a = 0. 

On making in (A.l) the change of variables 

w = K v = - (A2) 

y y 

we get 

w = w 2 + k, v + kv = 2[3y 3 . (A3) 

First, we consider the case k = —a 2 < 0. Then the general solution of the first equation 
from (A.3) is given by one of the formulae 

w = —acotha(t + Ci), w = — atanha(t + Ci), w = ±a, C\ E R 1 , 

whence 

y = C 2 sinh _1 a(t + Ci), y = C 2 cosh" 1 a(t + d), y = C 2 exp(±at) , C 2 G R 1 . (AA) 

The second equation of system (A.3) is a linear inhomogeneous ODE. Its general 
solution after being substituted into (A. 2) yields the following expression for z(t): 

(C 3 cosh at + C 4 sinh at) sinh" 1 a(t + C\) H ^- sinh -2 ait + Ci), 

/3C 4 

(C 3 cosh at + C 4 sinh at) cosh" 1 ait + CA + cosh" 2 ait + CA, ( A5) 

a 1 

(C 3 cosh at + C 4 sinh at) exp(±at) + — — |- exp(±4at), 63,64 C i? 1 . 

The case fc = b 2 > is treated in an analogous way, the general solution of (A.l) being 
given by the formulae 

y = D 2 sin" 1 6(t + £>!), 

3D 4 ( A6 ) 
z = (D 3 cos bt + D 4 sin bt) sin -1 6(t + D x ) + -— - sin" 2 bit + DA, 

b l 

where D ± , D 2 , D 3 , D 4 are arbitrary real constants. 
The case 2. a^O, (3 = 0. 
On making in (A.l) the change of variables 
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we have 



z 

y = exp w, v — - 

y 



id — w 2 = k + aexp4u>, v + kv = 0. (A7) 
The first ODE from (A. 7) is reduced to the first-order linear ODE 

1 dp{w) 



2 dw 

1/2 



— p{w) = k + a exp Aw 



by the substitution w = (^(w)^ , whence 

p(w) = a exp Aw + ^exp2w — k, 7 E R 1 . 

1 /2 

Equation w = (p(w)^j has a singular solution w = C = const such that p(C) = 0. 
If w 7^ 0, then integrating the equation w = p(w) and returning to the initial variable y 
we get 

vit) 

( dT -t\r 

J r(ar 4 + 7 r 2 - k) 1 / 2 1- 

Taking the integral in the left-hand side of the above equality we obtain the general 
solution of the first ODE horn. (A.l). It is given by the following formulae: 
under k = —a 2 < 

y = C 2 (a + sinh 2a(t + dj) ^ , 

y = C 2 (a + cosh2a{t + d))^ 2 , (A8) 



y = C 2 (a + exp(±2at 



-1/2 



under k = b 2 > 



y = D 2 ( a + S in 2b(t + D 1 )) 1/2 . (A9) 
Here C±, C 2 , D±, D 2 are arbitrary real constants. 

Integrating the second ODE from (A. 7) and returning to the initial variable z we have 
under k = —a 2 < 

z = y(t)(C 3 cosh at + C 4 sinh at) (A 10) 

under fc = b 2 > 

z = 2/(t)(D 3 cosfet + D4sin6t), (All) 
where C 3 , C 4 , D 3 , D 4 are arbitrary real constants. 
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Thus, we have constructed the general solution of the system of nonlinear ODEs (A.l) 
which is given by the formulae (A.5)-(A.ll). 

Appendix 2. Separation of Eq.(|T|) into three second-order ODEs. 

Suppose that there exists an Ansatz 

u = Q(t, x)ip (u (t, x^ip^u^t, x)^if 2 {uJ2(t, x)^j (A. 12) 

which separates the Schrodinger equation into three second-order ODEs 



d?ip 
dujQ 


= U (u ,ip , 


dipo 

dido 


d 2 (pi 


= Ux(ui,(pi, 


d(fi 


duj\ 


duii 


d 2 ip 2 


= U 2 (uJ 2 ,ip 2 , 


dip 2 


dw 2 


duj 2 



; Al, a 2 ) 



(A.13) 



according to the Definition 1. 

Substituting the Ansatz (A. 12) into Eq.(|I]) and excluding the second derivatives -j-tt, 

______ ^ 

fi = 0, 2 according to Eqs.(A.13) we get 



+2Q Xa uj 0xa ip (pi(p 2 + 2Q Xa u 1Xa (p (p 1 ip 2 + 2Q Xa uj 2Xa (p ipi(p 2 + Q((Au )<p (pi(p 2 

+ (Aa;i)v9 ^lV 3 2 + (AW 2 )V>0</>1<P2 + ^Ox^OxJJq^I^ + ^\x a ^\x a <-P^J\^2 
+^2x a ^2x a ( P0 ( PlU 2 + 2uJ 0Xa UJi Xa ip ipiip 2 + 2u OXa UJ2 Xa (po(pi<p2 + ^lx a ^2x a 'PQ < Pl < P2) 

= VQip fiip2- 

Splitting the above equality with respect to 0oy^i> <fio<p2, ^1^2 we obtain the equalities: 

uox a Ui Xa = 0, uj 0Xa u 2Xa = 0, uj 1Xa u 2Xa = 0. (A. 14) 

Since the functions oj^, /i = 0,2 are real-valued, equalities (A. 14) mean that there 
are three real two-component vectors which are mutually orthogonal. This is possible 
only if one of them is a null-vector. Without loss of generality we may suppose that 
(u 0xi , u 0x2 ) = (0,0), whence u = f(t) ~ t. 

Consequently, Ansatz (A. 12) necessarily takes the form (|5]). But Ansatz (|5p can not 
separate Eq. (fj) into three second-order ODEs, since it contains no second-order derivative 
with respect to t. 

Thus, we have proved that the Schrodinger equation ([]]) is not separable into three 
second-order ODEs. 
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